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Q ■ Abstract. A new approach in the quantum theory of few-electron nanoelectronic devices 

- the S-matrix approach - is presented in a simple example: a single-electron transistor 
consisting of a single-level quantum dot connected with two metallic leads through the 
corresponding potential barriers. The electron transport through the quantum dot due 
to the electron tunneling between the dot and the leads is studied. The strong Coulomb 
' repulsion between the electrons in the dot is taken into account exactly, while the tunneling 

r**"* ■ between the dot and the leads, considered as a small perturbation, is studied by means of the 

perturbation theory. For summing up the infinite perturbation theory series we apply the 
Green function technique and the Heisenberg equation of motion of the electron annihilation 
and creation operators. The matrix elements of the transition processes include both the 
. direct and crossing terms, so that there is no need to use the non-crossing approximation 

(NCA). The explicit expression of the electron transport current is derived. 
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Ph! I. INTRODUCTION 

The transport of electrons through the single-electron transistor (SET) consisting 
of a quantum dot and two (metallic or ferromagnetic) leads under some bias voltage 
between two leads was studied in many theoretical and experimental works. Following 
Meir, Wingreen and Lee [1] for the study of the electron transport many authors [2- 11] 
used the Landauer formula [12] expressing the conductance in terms of the Green function 
of the electron in the dot, while Izumida et al. [13- 15] applied the Kubo linear response 
theory [16] to express the electrical current (or the conductivity tensor) in terms of the 
Green functions. The calculation of the Green functions has been done by means of 
different approximate methods: the truncation of the system of equations for many-point 
Green functions derived from the Heisenberg equation of motion [1,3], the perturbation 
theory with respect to the strong Coulomb repulsion potential [2,8,9], the non-crossing 
approximation (NCA) [7], the numerical renormalization group method [10,13,15] and 
the Quantum Monte Carlo technique [14]. 

In a recent work [17] for the study of the electron transport through a single-level 
quantum dot one of the authors has proposed another method based on the use of the 
S-matrix in the perturbation theory with respect to the Hamiltonian of the tunneling 
between the dot and the leads considered as the interaction Hamiltonian Hint, the strong 
Coulomb repulsion potential between two electrons in the dot being included into the part 
Hq with exactly determined eigenvalues and eigenvectors of the total Hamiltonian 



H — Hq + Hi n t ■ 



(1) 
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As in earlier works [1 - 11, 13- 15] we consider the single-electron transistor consisting 
of a single-level quantum dot connected with two metallic leads "a" and "b" through two 
potential barriers. The dot and the leads are assumed to be kept in the quasi-equilibrium 
states at a fixed temperature and with given chemical potentials being under control. We 
use the unit system with h = c = 1 . 

Denote c CT and c+ the annihilation and creation operators of the electron with the 
spin projection a =j, I at the single energy level E in the quantum dot, a a (k), b a (k) and 
a+(k), b+(k) those of electrons with the spin projection a and the momentum k in the 
leads "a" and "b", resp. Following earlier works [1 — 11, 13 — 15] we assume that 

Ho = Ej24c a + Un^ ni + ]T ]T { /T„(fc) </,; ifc) a a (k) + E,,(k) b„ (*) Mfc) } , (2) 

cr k cr 

n a = c+ c CT , (3) 



H int = E E { V *W a - W C - + Va ^* C ° °*( fc ) + ^ Ca + V b (kY C+ b a (fc)} . 

k cr 

(4) 

In the previous report [17] the electron current from one lead to another one in the 
presence of some bias voltage between two leads was calculated in the second order of the 
perturbation theory with respect to the interaction Hamiltonian (4). In this work the sums 
of the series of high order terms of the perturbation theory in the ladder approximation 
are calculated by means of the Green function technique with the use of the Heisenberg 
equation of motion for the electron annihilation and creation operators. As we shall see, 
the matrix elements of the S-matrix contain both the "direct" and "crossing" terms, so 
that in the framework of the present S-matrix approach there is no need to assume the 
non-crossing approximation. The "crossing" terms describe the co-tunneling processes. In 
many cases their contributions are significant. 

II. PERTUBATION THEORY IN THE SECOND ORDER 

In order to apply the perturbation theory we introduce the electron annihilation and 
creation operators as well as the interaction Hamiltonian in the interaction representation 

c+(t) = e iHot 4 e~ iHt> \ c a (t) = e iHot c a e- iHo \ 

a+(Jb, t) = e iHot a+(Jb) e~ iH °\ a a {k, t) = e iHot a a (k) e~ iHot , 

b+(k, t) = e iHot b+(k) e- iHot , b a (k, t) = e lHo %{k) e~ iHot , (5) 
H mt (t) = e iHot H int e- iHot . 
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The S-matrix in the perturbation theory with respect to the interaction Hamiltonian (4) 
equals 

oo oo oo 

S = rjexp (-i J H int (t)dtj^ = E^T~ J dtl - J dt n T{H int (t l )...H mt (t n )}. 

-oo n=0 -oo -oo 

(6) 

S=l-iR, (7) 



We set 
and write 



R = J2R {n \ (8) 

i=i 

where denotes the contribution of all n-th order terms. In the first order of the 
perturbation theory we have 

oo 

= EE / dt { V a (k)a+(k,t)c a (t) + V a (k)*c+(t)a a (k,t) + (9) 

<T k — OO 

+V b (k) b+(k, t) c a (t) + V b (k)*c+(t) b a (k, t)} , 

while in the second order 

OO oo 

* (2) =-;;eeee/^ j dt2 

CTl CT2 k! k 2 _oo _oo 

T { [ V a (fci ) o+ (fci , 1 1 ) c CTl (t x ) + V a (fci ) *c+ (ti ) a CTl (fci , 1 1 ) 

+V 6 (fci) 6+ (fci, ti) c CT1 (tx) + H(*i)*c+ (ti) (*!,*!)] (10) 

[ K(fc 2 ) 0+ (fc 2 , t 2 ) C CT2 (t 2 ) + K(fc 2 )*C+ (t 2 ) CL(j 2 (* 2 , t 2 ) 

+Vt(fc 2 ) 6+ (*2, *2) c CT2 (t 2 ) + H(*2)*C+ (t 2 ) 6 CT2 (fc 2 , i 2 )] } . 

If the initial and final states |«) and |/) in some matrix element (/ | R \ i) of the 
scattering operator i? contain one and the same subsystem of n particles a a (k) and m 
particles b a (k), (n,m=0,l,2,. . . ) which do not participate in the interaction process, then 
this matrix element can be written in the form 

</|i?K> = ^(-l) p <a < (fc / 1 ) \a tT1 (k 1 ))---(a < (k' n ) \a an (k n ))x (11) 
p 

)){f\R\i), 

where { i } and | / | are the subsystems complementary to those of the non-interacting 
particles a ai (ki), ... a CTn (fc n ), b an+1 (k n+1 ), ... b a 

n+m (fcn+m) m the initial and final states, 
resp., P denotes the antisymmetrization with respect to the identical particles of each 
kind a a (k) and b a (k). The matrix elements in this factorized and antisymmetrized form 
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are called reducible ones contrary to the irreducible matrix elements which cannot be 
expressed in the above-mentioned form (11). In the study of the observable physical 
quantities of the system it is necessary and sufficient to derive the explicit expressions of 
the irreducible matrix elements of the scattering operator R. 

Consider now in detail the transport of the electron through the quantum dot - the 
transition of the electron from one lead into another one (a a (k) — > b a i{k') } for example) 
via the intermediate states of the electron in the quantum dot (virtual particles c CT ) by 
means of the tunneling of the electrons between the leads and the dot. In the second order 
of the perturbation theory we have following irreducible matrix elements: 

(b a ,(k>)\RW\a a (k)) = 

OO OO 

- iv a (k)*v b (k') f dh f dt2 ^ b {k)t x -E a {k)t 2 ] { q |r{ c<T ,(*i)c+(t 2 )}|o>, [LZ) 



■ (T2 
CO 



,b a[ (k')\RW\a ai (k)c 

■■V a (k)*V b (k') f dh f ^ 2e ^(fc')*i-^(fc)* 2 ]/ C(T , TfoWc+famco 



(13) 



(c i c 1 b af (k')\R^\a a (k)c^c l ) = 

oo oo (*\ A\ 

- i V a (k)* V b {k') f dh f dt 2 e W* )ti-E a (k)t 2 ] (qCT \ T {c a ,(h) c+(h)} |c T q ). ^> 

— oo — oo 

The problem is now reduced to the calculation of the matrix elements of the chronological 
product of two operators c a i(t') and c+(t) between different electronic states of the quan- 
tum dot: the empty state |0), the singly occupied states \c a ) and the doubly occupied one 

|c T q)- 

The matrix element 

G^ a (t'-t) = -i(0 \T{c a ,(t')ct(t)}\ 0) (15) 

is the usual Green function of the electron in the quantum dot, while two other matrix 
elements 

(16) 
and 

G { X(t' -t) = -i <qc T \T{c a ,(t') c+(t)}\ CT q> (17) 

are only similar to the usual Green function and might be called the generalized Green 
functions. They can be exactly calculated without any approximation with respect to the 
large Coulomb repulsion potential U. 

The chronological product of two operators c a i{t') and c+(t) has the explicit form 

T{(v(0 c+(t)} = 6(t' - t) <v(0 c+(t) - 8{t - c+(t) <v(0 . (18) 



,0) 



GT , it' 



t) 



T{c a ,(l/)c+(t)} 
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The first term in the r. h. s. of the relation (18) describes the creation of an electron 
in the quantum dot at the moment t (due to the tunneling of this electron from the lead 
"a") and then the subsequent annihilation of this electron {a = a') or another electron 
(a ^ a') in the quantum dot at a later moment t' > t (due to its tunneling to the lead 
"6"), while the second term describes the annihilation of an electron in the quantum dot 
at the moment t' (due to its tunneling to the lead "6") and the subsequent creation of 
another electron in the quantum dot at a later moment t' > t (due to the tunneling of this 
electron from the lead "a"). The first term is called the direct one, while the second one 
is called the crossing term. They are represented by the diagrams in Fig. a and Fig. lb. 

t'> t v t> t' 

b 




Fig. la 



Fig. lb 



The ability to take into account the contributions of the crossing terms is an advan- 
tage of the present S-matrix method. If the initial state is the empty one |0), then there is 
no electron in the quantum dot to be annihilated and the crossing term does not appear: 



<0 \T{c a ,(t') 4(t)}\ 0} = 9(t> - t) e- iE V-V5 a(T , 



(19) 



On the contrary, if the initial state is the doubly occupied one |c-|-cj_), then it is impossible 
to add one more electron into the quantum dot and only the crossing term gives the 
contribution: 



<qc T \T{c a ,{t') 4{t)}\ c T q) = - 6(t - t') e'^+V ^S aa > 



(20) 



In the case of the singly occupied initial states \c a2 ) both the direct and crossing terms 
can give non-vanishing contributions depending on the spin configuration of the electron 
system: 



CTi <J2^—cr' 1 a' 2 



(c. a , \T{c^) c+ (t)}\ c a2 ) = 9(t' - t) e-« E +^ 

(21) 

Introduce the Fourier transformations of the usual and generalized Green functions 



1 



(t' -t) = — 
' 2vr 



(22) 
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G#, (t'-t) = -!- [due-**? (u), (23) 

G ( X(t' - t) = ±- J due-^-^G%(u). (24) 

From the expressions (|19[) - (21) of the Green functions it follows that their Fourier trans- 
forms equal 

- df^ • < 25 > 

4<7i<Ti CT 2 ^ ~ w -E-U + io u-E-io ' 

= ^ . • (27) 
uj — h — U — 10 

It is easy to express the matrix elements (13) - (14) in terms of these Fourier trans- 
forms. Finally we obtain following relevant non-vanishing irreducible matrix elements 
contributing to the transitions of a particle a a (k) into a particle b a /(k') via the intermedi- 
ate virtual particle c a » in the second order of the perturbation theory: 

<6 T (A:') |^ 2 >|a T (fc)> = (^(A:') | ^( 2 ) | a x (fc)> = 

= 2vr 5[E a (k) - E b (k')] V a (k)*V b (k') • , . , (28) 

< Ci 6 T (fc')|^|a T (fc) Ci > = (cMk')\R {2) \^k)c,) tj + l ° 

= 2n 5[E a (k) - E b {k')\ v a (kyv b (k') ■ — — L— . , , ( 29 ) 

<c T 6 T (fc') \R^\a,(k) c T > = (c lbl (k>) iR^a^fL U + % ° 

- 2vr 6[E a (k) - E b (k')) V a {k)*V b {k') ' (30) 



[cMk')\R^a,(k)c l ) = (cMk')\R^\a l (k)c,)= Ea{k) E l ° 
= 2tt 5 [E a (k) - E b (k')} V a (kyV b (k') 



1 1 

+ 



(31) 



E a (k) - E - U + to E a (k)-E-io 

(c lC Mk')\R i2) \a r (k)c m ) = (c^b^lR^la^k)^) = 

= 2vr 5[E a (k) - E b (k')} V a (k)*V b (k') ■ 1 , • (32) 

h a {K) — h — U — to 

Note that the matrix element ()31|) describes the spin-flip scattering process inside the dot 
while other matrix elements are those of the non-spin- flip ones. 

Using the expressions (28) - (32) of the matrix elements of the scattering operator 
it is easy to calculate the transition rate. If the matrix element of the transition from an 
initial pure state \i) to a final one |/) has the form 

(/ | R | i) = 2vr S[Ei - E f ] M^ f , (33) 
then the transition rate in an unit time interval equals 

W i ^f = 2irS[E i -E f \ \M^ f \ 2 . 
If the initial state is a mixed equilibrium state / containing the pure state \i) with a weight 
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Wi , then the transition rate from this mixed state I to the final states |/) belonging to 
some class F equals 



iei feF 



ta. 



(34) 



Let us apply this formula to the transitions between the initial and final states 
presented in the matrix elements (28)- (32). The statistical weights of these pure states 
in the initial mixed equilibrium state are given in the Table I with 



Z — Z a Z\) Z c , 

Z a = [1 + e-WO] 2 , 

Z 6 = [l + e-W*)] 2 , 

Z c = l + 2e~P E ' + e-^ 2E '+ u \ 



Table I 



(35) 



State 


Weight 


\a c (k)) 


e -PE> a {k)/ z 


M fc )°i( fc )) 


e -(12E' a (k) /Z 


\a a {k)b a ,{k')) 


e ~/3[E' a (k)+E' b (k')y z 


|a T (*)M*') «!(*)> 


e -p[2E> a (k)+E> b (k')]/ Z 


\a a (k) <v) 


e ~f3{E> a (k)+E']/ Z 


|a T (fc)<v a l( k )) 


e -f3[2E> a (k)+E>]/ Z 


\a a (k)c a ,> b a/ (k')) 


e -p[E> a (k)+E> b (k')+E>}/ z 


\ai(k) c„ii b a r(k') ai(k)) 


e -(3[2E> a (k)+E b (k')+E'}/ Z 


\a a (k)ci q) 


e -/3[E' a (k)+2E'+U] j Z 


|a T (A;)c T cj. a^Jfe)) 


e -p[2E' a (k)+2E>+U]/ Z 


|a CT (A;)c T q 6 CT ' (A/)) 


e -f3[E' a (k)+E b (k')+2E'+U] j Z 


M fc ) c t cj. 6 (7 /(fc , )a|(ib)) 


e -f3[2E' a (k)+E' b (k')+2E'+U] j Z 
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where 



E' a (k) = E a (k)- f i a , 
E' b (k) = E b (k) - fi b , 
E' = E-fi c , 



(36) 



fi a , [j>b and fj, c being the chemical potentials in the leads "a", "b" and in the quantum dot, 
resp. Substituting the expressions of the matrix elements (28) - (32) and the statistical 
weights in the Table I to the r. h. s. of the formula (34) we obtain following expression of 
the probability of the transition of electrons from the lead "a" to the lead "b" through its 
intermediate states in the quantum dot: 



oo 



1 



— oo 



(37) 



W(u) 



/N |2 



J2^S[u-E a (k)} \V a (k)\ 2 Y J ^^^- E b{k')\ \V h {k')\ 

k k' 

1 2 / 



1 



+ 



lo - E + io 
1 



-pE' 



lo - E -U + io 

2\ 



+ 



lo — E — U + io lo — E — io 



+ e -/3(2E'+U) 



lo — E — io 
1 



+ (38) 



lo — E — U — io 



Here f a (co) and /b(w) are the Fermi distribution functions of the electrons in the leads "a" 
and "b", resp., 



I -)_ e -/3(uj-n b ) 



■ /.(w) = 



e -P(w-Ha) 
I _|_ e -P{u-Va) 



(39) 



For the reverse transition of electrons from the lead "b" to the lead "a" we have a similar 
expression of the transition rate 



1 



7T 



dwAM [i-/„H] w(lo) 



(40) 



If /U a > /ib, then W a ->b > Wb^ a and the effective transition of electrons from the 
lead "a" to the lead "b" generates the electrical current 



oo 

J e = e[W a ^ b - W b ^ a ] =e^ J dio[f a (Lo)-f b (Lo)] W(uo) 



(41) 



from the lead "b" to the lead "a" , where e is the absolute value of the electron charge. 

In the expression (38) there are the contributions of both direct and crossing terms 
of the chronological product (21). If we neglect the crossing terms, then we have the non- 
crossing approximation (NCA). In this approximation W(lo) is replaced by the following 



Quantum Theory of Electron Transport Through Single-Level Quantum Dot 



9 



expression 



W NC a(u) =Y,^5[u-E a {k)\ \V a (k)\ 2 Z^Hu-Ehik')} \V b (hf)r x 
fe fc' 



i 



(jj—E+io 



u-E-U+io 



(42) 



Obviously, if e is not very small, then the contribution of the crossing terms might be 
comparable to that of the direct ones. 

III. HIGH ORDER APPROXIMATIONS OF 
THE PERTURBATION THEORY. GREEN FUNCTION TECHNIQUE 

The expression in the r. h. s. of the formula (38) exhibits the sharp resonances at 
the values of oj near the energy levels of the electron in the dot u ~ E and u ~ E + U. 
However, at these values of u the broadening of the energy levels due to the tunneling of 
the electrons from (to) the dot to (from) the leads must be taken into account. This can be 
done by summing up the contributions of the high order terms in the ladder approximation. 
The most convenient way to do that is the application of the Green function technique. 

As an example consider again the transport of the electrons from the lead "a" to 
the lead "b" due to the single-electron tunneling through the quantum dot. Instead of 
the annihilation and creation operators a a {k), b a (k) and a+(&), b^(k) for the electrons 
in the leads in the Schrodinger representation the S-matrix (6) is expressed in terms 
of the corresponding operators a a (k,t), b a (k,t) and a+(fc,t), b^(k,t) in the interaction 
representation defined according to the formula (J5J) with the expression (2) of the operator 
Hq . It is easy to verify that 



(43) 



a+(Jb,t) = e iEa{k)t at(k) , a a (k,t) = e'^^a^k) , 
b+ (Jfe', f) = e ^( fe ')*' 6+ (&') , b(ji (k', t>) = e-^^'V (k') , 

From the expression (6) of the 5-matrix 

oo 

5 = r|exp (-i J H int (t)dt\\ 

—oo 

with Hi n t(t) being determined by the formulae (4), © and (43) 

Hintit) = EE{ V a (k) e^-W'o+C*) c a (t) + V a (k)* e-^^ctit) a^k) 

k a (44) 

+V b (k) e lE ^ k Ht(k) c a (t) + V b (k)* e-^^c+it) b a (k)} 

it follows that 

(MAOISM*)) = 

- iV a (k)*V b (k') J df J dt^mt'-Eai^t]^ \ T {c a ,(t')c+(t)S}\0), (45) 
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\Ca' 2 K[(k') \R\a ai {k)c a2 ) = 

-iv a {kyv b (k') J dt> J ^e 4 ^( fc ')*'~^( fc )*](c 4 |r{ C(T ,(Oc+(t)s}| C(72 >, (46) 



-oo — oo 



(cjcffv^) \R\ a a (k) C|C|) = 

- % V a (k)* V b {k') J dt' J dte i [^(*') t '-^( fc ) t ]<c i c T \T{ Cu ,{t') c+(t)S} |c T q >. (47) 

— oo — oo 



The calculation of these matrix elements (45) - (47) requires the determination of following 
quantities 

G^ c {t' -t) = -i(0 \T { cAt') 4(t)S } | 0) , (48) 
G^^'-t) = -^|r{c fri (t')c+(t)S'}|c CT2 > ) (49) 

G { X CC (t' ~t) = -i(c^\T { <v(0 c+(t) S } | c T q> (50) 

which might be called the (generalized) Green functions of the electrons in the dot in 
the presence of the tunneling transitions to and from the leads. The expressions of these 
(generalized) Green functions may be derived by means of the Green functions technique. 

In order to apply the Green function technique we must work in the Heisenberg 
representation. The transformation from the interaction representation to the Heisenberg 
one is performed with the use of the operator S(t,to) satisfying the Schrodinger equation 
with the interaction Hamiltonian Hi n t(t), 

i^^- = H int (t)S(t,t ), (51) 

and the initial condition 

S(t ,t ) = l. (52) 

The solution of the equation (|51|) with the initial condition (|52j) can be represented in the 
form similar to the formula (6) for the S-matrix: 

t> 

(53) 



S(t',t ) = rjexp (-i J H int {t)dt\ | 
to 



It is obvious that 

S = S(oo, -oo) . (54) 
From the expression (|53|) it follows the group property of the operators S(t,to): 



S(t,t ) = S(t,t')S(t',t ) 



(55) 
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From this property and the initial condition 1)52(1 we obtain the relation 

S(t ,t) = S(t,t )- 1 . (56) 
Moreover, the operator 1)53|) is unitary 

S(Mo) + = S(Mo)- 1 =S(to,t). (57) 

Now we derive the explicit expression of the chronological product T { c a > (t') c+ (t) S } 
in the matrix elements (|48|) - (|5U() . In the case t' > t we write 

S = 5(oo, — oo) = S(oo, t) S(t', t) S(t, — oo) 

and have 

r{tv(0 ^} = s(oo, <v (0 s(t', *) c a (*) s(t, -oo) 

= 5(oo, -oo) S(-oo, t') c a , (f) S(tf , -oo) 5(-oo, t)c+(i) S(t, -oo) . 



Introducing new operators 



C a ,(t') = S{-oc,t')c a <{t') S(t', -oo), 
C+(t) = S(-oo,t)c+(t)S(t,-oo) 



(58) 



we rewrite 

T{c a ,(t')ct(t)S} = SC a ,(t>)C+(t), t'>t. (59) 
Similarly, in the case t > t' we have 

T{c af (i/)c+(t)S} = -SC+(t)C a r(t'), t>t'. (60) 

Combining both expressions (|59|) and ()6(J|) finally we obtain the formula 

T{ C<J ,(t')c+(t)S} = ST{C a ,(t')C+(t)} . (61) 

Now the (generalized) Green functions (|48[) - (|50[) are expressed in terms of the new oper- 
ators <|55)l : 

- *) = "*<0 |Sr{C ff ,(0 C+(i) } |0> , (62) 

G^ ai Jt' -t) = -i(c a ,\ST{ C a ,(t') C+(t)}\ C<72 ) , (63) 

G^{t' -t) = -i{c lC] \ST{C a ,{t') C+(t)}|c TCi > . (64) 

In order to derive the differential equations for the (generalized) Green functions 
(62) - (64) we must use the Heisenberg equations of motion for the new operators (|58|). 
Denoting 

H(t) = S(-oo, t) [H + H int (t)] S(t, -oo) (65) 
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the total Hamiltonian in the new representation we obtain the Heisenberg equations of 
motion with the total Hamiltonian H(i): 

i ^l = [ C a{ t),H(t)), i^^ = [C+(t),H(t)]. (66) 

This means that C a (t) and C£ (t) are the annihilation and creation operators in the Heisen- 
berg representation. 

Together with the operators C a (t) and C^(t) determined by the formula in- 
troduce also the annihilation and creation operators for the electrons in the leads in the 
Heisenberg representation 

A a (k, t) = S(-oo, t) a a (k, t) S(t, -oo) , 
A+(k, t) = S(-oo, t) a+(k, t) S(t, -oo) , 

B a (k,t) = S(-oo, t) b a {k, t) S(t, -oo) , { ' 

B+(k, t) = S(-oo, t) b+(k, t) S(t, -oo) . 

They satisfy also the Heisenberg equations of motion with the total Hamiltonian H(i): 
JA^t) = t)R {t)] l d -^A = [A +(k, t),H (t)] , 

i d_Bj±t± = [BAM)H{t)h = [B+(k,t),H(t)]. 

The total Hamiltonian in the Heisenberg representation is expressed directly in 
terms of the new operators (|58j) and (|67|) : 

H(t) = E^C+(t)C a (t) + UN 1 (t)N l (t)+ 

+ EE {E a (k)A+(k, t) A a {k, t) + E b (k)B+(k, t)B a (k, t)} + 

k * (69) 
+ EE(^„ (k)A+(k, t) C a {t) + V a (k)*C+(t) A a (k, t)+ 

+V b (k)B+(k, t) C a (t) + V b (k)*C+(t) B a (k, t )} , 
N a (t) = C+(t) C a {t) . 

The operators (|58|) and l|67|) satisfy the canonical equal-time anticommutation relations 

{C a (t), C+(t)} = 5 aa ,, 

{A a (k,t) , A+{k',t)} = 6 ( , a ,5 w , (70) 
{B (k,t) , B+(k',t)} = 5 aa ,5 kk ,, 

other equal-time anticommutators between these operators being equal to zero. 

From the Heisenberg equations of motion (|66[) and (j68j) . the expression (|69|) of the 
total Hamiltonian, the equal-time canonical anticommutation relations (|7U|) and other ones 
with the vanishing r. h. s. we derive following differential equations for the operators (|58|) 
and (JHZJ): 



. dC a (t) 
1 dt 



EC a (t) + UN. a (t) C a (t) + i Va(k)*A a (k, t) + V b (k)*B a (k, t)} , (71) 
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i dA °^ t] = E a (k)A a (k, t) + V a (k) C a (t) , (72) 

i dB ^ t] = E b (k)B a (k, t) + V b {k) C a {t) . (73) 
at 

Using these differential equations and the canonical equal-time anticommutation relations 
we can derive the systems of equations for the (generalized) Green functions. 

IV. CALCULATION OF THE (GENERALIZED) GREEN FUNCTIONS 

For the convenience in the determination of the real and imaginary parts of the 
Fourier transforms of the (generalized) Green functions (|62|)— (64) we divide each of them 
into two functions according to the definition of the chronological product 

T { Ci,t>) C+(t) } = e(t> - t) C a ,(t')C+(t) - 9(t - o C+(t) C a ,(t>) . 

We have 

G ( T(t' ~t) = G ( l Cci+) (f -t) + «<-> (tf - t) , (74) 

GST™ (f - *) = G®;£1 [f -t) + G^J<;1 {f - t) , (75) 

G<?«(f - t) = G^ c4+) (f -t) + G^ CC{ - ] (f " t) (76) 

with 

G (o)cc(+) {t , _ t) = _ i9{t > _ t) ( Q \ SCal (t') C+(t) | 0) , (77) 

Gj;'"^ ~t)= i0{t - t')(0\SC+(t) 0^(0 I 0} , (78) 

G^l(t' -t) = -i0{t> - t)(c a ,\SC a[ (t')C+(t)\c a2 ) , (79) 

G^^Jt' -t) = i9(t - 0(^1^(00^(01^) , (80) 
G (2)cc(+) (t , _ t) = _ iQ{t , _ t) ^ Ci \ SCAtl) C+(0|qc T > , (81) 
G P)«H (t / _ t) = i6{t _ t'^c^SC^CA^M ■ (82) 

Using the equation (71) for 0^/(0 we can derive the differential equations for the 
functions G^)^ c ^ (t' — t) , G,,^ (t' — t) and O^^f (f — t) which contain many other 

O 2 j 1 2 

(generalized) Green functions. The derived differential equations do not form closed finite 
systems of equations, but belong to some infinite systems. In order to find their approx- 
imate solutions we might truncate them at a corresponding order and obtain the closed 
systems of approximate differential equations which can be solved exactly. Let us now 
truncate the infinite systems of equations by neglecting the 4-point (generalized) Green 
functions multiplied with small functions V a (k), Vb(k), V a (k)*, Vb(k)* and all n-point 
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equations containing following time-independent matrix elements (constants) 

(o\sc (7 ,(t)c+(t)\o) = g ^ cc{+ \ 

(0\SC+(t)C a ,(t)\0)=g^-\ 

(cvj5^ i (t)C^(*)|^)=^ + Ji i , 

(ca il \SC+(t)C oi (t)\c n )=g^^ 

(c i0t \SC < ,(t)C+{t)\c r c i ) = g { X CC(+ \ 

(c^\SC+(t)C al (t)\c^c i )=g^ CC{ ~\ 
(0\SA a ,(k,t)C+(t)\0)=g^ + \k), 

(0\SC+(t)A a ,(k,t)\0) = gf): c{ -\k) , 
(c a , \S A a ,(k,t)C+(t)\ CtT2 ) = g® a £i(k) , 
(c a ,\SC+(t)A a ,(k,t)\c a2 ) = - 
< Ci c T \SA a ,(k,t)C+(t)\ c TCi > = g { Z ac{+) (k), 
(qc T \SC+(t)A al (k,t)\ c lH ) = g { X ac{ -\k) , 
(0\SB a ,(k,t)C+(t)\0)=g^ +) (k), 
(0\SC+(t)B a ,(k,t)\0) = g$W-\k) , 
(c a ,\S B a ,(k,t)C+(t)\c a ,) = 9 ® b ff 2 (k), 

(c a ,\SC+(t)B a ,(k,t)\c^ = 9 ® b ; { J 2 (k) , 
< Ci c T \SB al (k,t)C+(t)\c^) = g { X bci+ \k) , 
(qc T \SC+(t)B a ,(k,t)\c lCl ) =g^-\k), 
(0\SN_ a ,(t) C a ,(t)C+(t)\ 0) = /S cc(+) , 
(o\S N_ a , (t) Apt (k, t) C+(t) | 0) = ffl ac{+) (k) , 
(0\S N^, it) B a , (k, t) C+ (t) 1 0) = (k) , 

(0\SC+(t) N^(t)C a ,(t)\0) = /S cc( - } , 
(0 I S C+ (t) 4,, (k,t)\0) = /^° C(_) (fc) , 

<0| 5 C+ (t) N—a 1 (t) Bi,k, t)\0) = /<?>**-> (*) , 

<c CT , |5 AL^ (t) A a[ (k, t) C+ (t) | c CT2 > = /g^W (fc) , 
(<V 2 1 5 iV_ CT , (t ) B<T , (*, t) C + (t) | c CT2 > = (A;) , 
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fty^ (*) 



(c a/2 \SCX(t)N_ a/i (t)B a ,(k,t)\c^ 

< Ci c T ^JV^^C^^C+Ct) I c T c 4 ) = /^ c(+) , 
<qc T |5^(t)^(*,t)C+(t)|c TCi > = f$™ + \k), 
< Ci c T |5 iV_ CT , (t) B a , (k, t) C+(t) | CT q> = j^ fec(+) (fe) , 
<qc T |5C+(t) JV-^Ct)^^)! c T q> = f { X CC{ ~ ] , 
<qc T ^C+W^Ct) ^(fc,*)! c TCi > = /^^(fc) , 
< Ci c T | S C+ (t) AL ff , (t) Bp/ (k, t)|c T q) = /S 6c( - } (*) . 

Performing the Fourier transformation of the differential equations we derive a sys- 
tem of algebraic equations which can be solved easily and obtain following results: 

G^ cc{± \u) = {u>-E-Z(u;± io)}- 1 ■ (84) 



„(0)cc(±) 
9 a' a 



+E 

k 



uj — E a {k) ± io 



+ 



; - Eb(k) ± io 



+ U p(0)cc(±) i 
' a' a 



Q®CC{±) , , = r _ E _ s( )} 
(To O"'. (71 (TO V ' L V yJ 



-1 



m CC (±) 

"(Tj °1 CT 1 



ffla C (±) 



(k) 



+ 



uj - E a (k) ± io uj - Eb(k) ± io 



(85) 



G^\u) = {u>-E-i:{u±io)y l 



£ CC(±) +E 



,(2)6c(±). 



(86) 



F (0)cc(±) 
a' a 1 



{^-^-[/-^(wiio)}- 1 




(0)6c(±). 



+ 



a; - E a (k) ±io uj - E b (k) ± io 



(87) 



<J 2 o\a\<J2 >■ 



{a; -£-£/-£( a; iio)}" 1 
,JTJcc(±) 

» (T2(T^(Tl(T2 / / 



O r (±) (fc) ^.fflMi) (fc) 



+ 



w - £? a (A;) ± io w - i?fe(fe) ± io 
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where 



/ir (±) +£ 



,(2)ac(±) 
Ja'a 



(k) 



+ 



,(2)bc(±) 
Ja'cr 



(k) 



Scj ± io) 



cj - i?a(A) ± io to - Eb(k) ± io 

\v a {k)\ 2 \vm\ 2 



+ 



lo - E a {k) ±io lo - Eb(k) ± 



(90) 



is the self-energy part of the electron in the dot due to its tunneling from and to the 
leads. The self-energy part £(u;) in the denominators of the (generalized) Green functions 
describes the broadening of the energy levels of the electron in the dot. 

The constants (83) can be calculated in the perturbation theory with respect to the 
tunnelling Hamiltonian Hi n t. In the lowest (zero) order we have following values: 



Jcr'cr 
f [EJ CC (+) 



fT W (k) = fT {±) (k) = 0, 



0, 



/9r (±) (*) = sW +) (*) 



,(2)cc(+) 
J a' a 



0. 



A2)ac(±) 
J a' a 
(0)cc(+) 

if (±) (*) 

JD«(+) 

0-, 2 (T 1 <J1<J2 

JD«<-) 

JDac(±) »v 

g mce(+) = Q 

if {±) (^) = 



,(2)«(-) 
Vcr 

r(2)ftc(±) 



(k)=f^r m (k) 



= 0. 

= 0, 



(91) 



6(7' (7 



JO)cc(-) 
9 a' a 



^-(71(72^-^4(^(71 -^(7 2 ) 
^(71(72^(72(71 ' 

ffljc( ± ) (fc) 



5 (2)cc(-) = 



0. 



Substituting these values of the constants into the r. h. s. of the relations (84) - (89) we 
obtain 



(lo) = 0. 



p(0)cc(±) 



(72(7^(71(72 V 1 

F^T { - ] (u>) 

(7 2 (7^(71(72 v ' 

p(2)«(+) 



J (7 1 (72/ 



E — U — T,(u + io) 



0, 



(92) 

(93) 

(94) 
(95) 
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gT {+ \") = 
gT { -\") = 

(7 2 (7 1 <Tl<72 v ' 



uj — E — T,(uj + io) 

o, 

~ uj- E -U - £(w + io) 



uj — E — T,(uj — io) ' 



/V(2)cc(+) 

G (7'(7 



(w) = 0, 



5(7' (7 



E-U-E(u-io) 



(97) 
(98) 
(99) 

(100) 

(101) 
(102) 



Summing up the two parts of each (generalized) Green functions G^ cc (u;) , (-vl^^ (u>) 
and (5^^ :c (a;) finally we have 



<5(7'(7 



G? c , c ( w ) 

(7 2 (7^(71<T2 v ' 

i(2)cc 



uj — E — E(o; + io) 

^-(71(72 ^-CTi^ (^(71 ~ ^(7 2 ) 



£ - C/ - E(w + io) 



+ 



^<7 2 £(7 2 (7l 



E-Z(u 



IO) 



u) — E — U — T,(uj — io) 

VI. CONCLUSION 



(103) 
(104) 
(105) 



The expressions (103) - (105) of the (generalized) Green functions with the contri- 
butions of the high order approximations are similar to the expression (25) - (27) of these 
functions in the lowest order of the perturbation theory. The appearance of the self-energy 
part Ti(uS) in the denominators of the (generalized) Green functions is the consequence of 
the tunneling of electrons between the dot and the leads which takes place in high order 
approximations. Instead of the expression (38) of the function W(u) determining the 
transport current now we have 

W(lj) = ^2iv8{uj-E a (kj\ \V a {k)\ 2 ^2^5[uj-E b {k')} |W)| 2 - 

k k' 

2 



+ 




+ e 



-/3E' 



uj-E-U- £(u> + io) 



+ e 



uj — E — £(u — io) 

P(2E'+U) 



1 



uj — E — U — H(uj + io) uj — E — Yj(uj — io) 

2 "I 



uj — E — U — T,(uj — io) 



(106) 



In comparison with other methods the S-matrix approach presented in this work has 
several advantages. The strong Coulomb interaction between the electrons in the dot was 



Quantum Theory of Electron Transport Through Single-Level Quantum Dot 



19 



taken into account exactly without the use of the perturbation theory. All crossing terms 
describing the co-tunneling are included automatically so that there is no need to assume 
the non-crossing approximation (NCA). We have derived the expression of the electrical 
current generated by the single-electron tunneling between the dot and the leads. The 
role of the multi-electron co-tunneling will be studied in a subsequent work. Moreover, it 
is easy to extend the presented method to the study of the multi-level quantum dots as 
well as to other nanosystems. 

Since the electron spin projections are explicitly indicated in the expressions of the 
matrix elements, it is convenient to apply the presented method to the study of the spin- 
dependent transition (transport) processes in the spintronic materials and devices. By 
including the interactions of the electrons with the electromagnetic field it is straight- 
forward to extend the presented method to the study of the photon-assisted transition 
(transport) processes and develop the S-matrix approach in the theory of nanophotonic 
materials and devices. 
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